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We present the first observation of WW + WZ production in the channel with an identified lepton and two 
jets in 2.7 fb _1 of integrated luminosity collected with the CDF II detector in pp collisions at y/s = 1.96 TeV. 
The signal is separated from the large background using matrix element calculations and is observed with a 
significance of 5.4cr. The WW + W Z production cross section is measured to be f 7.7± 3.1 (stat)±2.4(sys) pb, in 
good agreement with standard model predictions. A complimentary measurement using a fit to the dijet mass 
is also presented. 



1. Introduction 

Measurements of heavy vector boson pairs (WW, 
WZ, and ZZ) are important tests of the electroweak 
sector of the Standard Model (SM) . Deviations of the 
production cross section from SM predictions could in- 
dicate the presence of new particles or interactions [l[ . 
Furthermore the topology of diboson events is similar 
to events where a Higgs boson is produced in asso- 
ciation with a W or a Z, so diboson searches and 
measurements are useful tests of analysis strategies 
employed in Higgs searches. 

WW and W Z production has been observed at the 
Tevatron in channels where both bosons decay lep- 
tonically 0. In semi-leptonic modes, where one 
boson decays to two quarks, large backgrounds make 
the signal extraction more challenging. CDF recently 
reported the first observation of such semi-leptonic 
decays in a channel with two jets and large missing 
transverse energy [ij • We present observation of semi- 
leptonic diboson decays in a channel with an identified 
lepton and two jets. The DO collaboration has previ- 
ously reported evidence of this signal in 1.07 Q. 

Discrimination of the signal processes from the large 
backgrounds is accomplished by use of matrix element 
probabilities. The implementation of this search tech- 
nique follows that used in the search for associated 
Higgs boson production Q . A second complimentary 
technique fits the invariant mass of the two-jet system 
(referred to as the dijet mass or Mjj). 

2. Experimental Apparatus 

The production is observed in 2.7 fb _1 of pp colli- 
sion data with yfs = 1.96 TeV collected by the CDF 
II detector, which is described in detail elsewhere Q- 
The aspects of the detector relevant to this analy- 
sis are described here briefly. The tracking system is 
composed of silicon microstrip detectors and open-cell 
drift chambers inside of a 1.4 T solenoid. Electromag- 
netic lead-scintillator and hadronic iron-scintillator 
sampling calorimeters segmented in a projective ge- 
ometry surround the tracking. The central calorime- 
ter covers \rj\ < 1.1 while the plug calorimeters extend 



the calorimetry into the region 1.1 < \r/\ < 3.6. Out- 
side of the calorimeter are muon detectors composed 
of scintillators and drift chambers. Cherenkov coun- 
ters around the beam pipe and in the plug calorimeter 
count the inelastic collisions per bunch crossing and 
provide the luminosity measurement. 

We use a cylindrical coordinate system with its 
origin in the center of the detector, where and 
4> are the polar and azimuthal angles, respectively, 
and pseudorapidity is r; = —In tan(0/2). The miss- 
ing E T 0t) is defined by $t = — ^2 % E^h il i = 
calorimeter tower number, where hi is a unit vector 
perpendicular to the beam axis and pointing to the 
i th calorimeter tower. $t is corrected for high-energy 
muons and jet energy corrections. We define the miss- 
ing transverse energy £r = \Bt I ■ The transverse mo- 
mentum pt is defined to be psmO. 



3. Event Selection 

The data samples used in this analysis are collected 
using trigger paths requiring a central electron (muon) 
with Et(pt) > 18 GeV as well as a trigger path re- 
quiring two jets and large $t- Offline we further 
select events by requiring an electron or muon can- 
didate with pt>20 GeV and exactly two jets with 
Et> 25 GeV. Jets are clustered using a fixed-cone 
algorithm with radius AR = y/ (Ar/) 2 + (Acf) 2 = 0.4 
and their energies are corrected for detector effects. 

Several event vetos are imposed to reduce back- 
ground levels and achieve good agreement between 
data and simulation. We reject events with any ad- 
ditional jets with E T >12 GeV. We reduce the back- 
grounds due to fully leptonic decays of it and diboson 
events by vetoing on additional charged leptons. We 
also reject events containing dilepton pairs consistent 
with decay from a Z boson, where the second lepton 
passes very loose identification cuts. Cosmic ray and 
photon conversion vetos are implemented. 

We impose additional requirements to reduce the 
level of background due to QCD multi-jet events. 
These events will enter our sample if a jet fakes an 
electron or muon and mismeasurement leads to large 



2 



Proceedings of the DPF-2009 Conference, Detroit, MI, July 27-31, 2009 



I£t- QCD multi-jet events are difficult to model, so 
we reduce their contribution as much as possible. In 
events with a muon, we require S! T >20 GeV and 
mxiW) >10 GeV, where mxiyV) is the transverse 
mass of the lepton-$r system. Jets are more likely to 
fake an electron than a muon, so we impose partic- 
ularly strict requirements in events with an electron, 
requiring 

$ T > 40 GeV and nvr(W) >70 GeV, as well 
as placing requirements on the $t significance and the 
angle between the -St and the second jet. 



4. Modeling and Backgrounds 

The dominant background to the diboson signal is 
VF+jets production where the W decays leptonically. 
Smaller but non-negligible backgrounds come from 
QCD multijet, Z+jet, tt, and single top production. 
The modeling for all processes except QCD multijet 
events is provided by event generators and a GEANT- 
based CDFII detector simulation. Pythia is used to 
model WW, WZ, and tt @ events. The W/Z+jets 
backgrounds are modeled using the fixed-order gen- 
erator Alpgen interfaced with the pythia parton 
showering framework Q. Madevent interfaced with 
Pythia is used for the single top background (loj . 
The QCD background is modeled using data events 
passing loosened lepton requirements. 

The normalization of WW, WZ, Z+jets, tt, and 
single top backgrounds is estimated using predicted or 
(in the case of Z+jets) measured cross sections and ef- 
ficiencies derived from simulation. The normalization 
of the QCD background is estimated by fitting the $t 
spectrum in data to the sum of all contributing pro- 
cesses, where the QCD and VF+jets normalizations 
float in the fit. While this fit provides an initial esti- 
mate of the V^+jets normalization, this normalization 
is a free parameter in the final fit used to extract the 
diboson cross section. The expected event yields are 
shown in Table Q] The observed yield agrees well with 
the total expected yield since the predicted yields are 
derived from a fit to the data. 

Table I Expected event yields for signal and background 
processes and total observed number of events. 



Process 


Event yield 


WW signal 


446 ± 17 


WZ signal 


79 ± 3 


W+jets 


10175 ± 152 


Z+jets 


584 ± 43 


QCD multijet 


283 ± 113 


tt 


131 ± 9 


single top 


110 ± 8 


Observed 


11812 



We expect 525 signal events out of 11812 total 
events. The dominant background is due to W+jet 
events. The QCD multi-jet background is small due 
to the stringent cuts imposed to reduce its size. 



5. Matrix Element Methodology 

A matrix element method is employed to separate 
signal and background events. The probability that 
an event was produced by a given process is deter- 
mined using the leading-order differential cross sec- 
tion for that process. For an event with measured 
quantities x, we integrate the appropriate differential 
cross section do~(y) over the partonic quantities y con- 
volved with the parton distribution functions (PDFs) , 
f(yi), /(j/2), and over transfer functions describing de- 
tector resolution effects, W(y,x): 

P(x) = da(y)dq 1 dq 2 f(y 1 )f(y 2 )W(y,x). (1) 

We use the CTEQ5L PDF parameterization [11]. The 
lepton momenta and jet angles are assumed to be mea- 
sured exactly. W(x, y) is a mapping of measured jet 
energy to partonic energy derived using the full de- 
tector simulation. The integration is performed over 
the energy of the partons and p^. The matrix ele- 
ment is calculated with tree-level diagrams from MAD- 
GRAPH [T3] • Event probability densities are calculated 
for the signal processes WW and WZ, as well as for 
the Wbb, Wcc, Wcj, Wjg, Wgg, and single top (i- 
channel and s-channel) backgrounds. 

The event probabilities are combined into an event 
probability discriminant: EPD = P s i gna l / (Paignal + 
Pbackground) , where P s i g nai = Pww + Pwz and 

Pbackground = P\Vbb + P\Vcc + Pw eg + P\Vgj + Pwgg + 

P st . We make templates of the EPD for all signal and 
background processes and ultimately extract the sig- 
nal using a fit of the observed EPD distribution to a 
sum of the signal and background templates. 

Figure Q] shows the dijet mass in bins of EPD. Most 
of the background events have low EPD, while events 
with EPD>0.25 are already very signal-like. The 
signal-to-background ratio improves with increasing 
EPD. 

Before comparing the observed EPD to the predic- 
tion, we validate the Monte Carlo modeling of the 
quantities that enter the matrix element calculation 
(the jet and lepton energies and angles). We divide 
events into three regions according to their Mjj: the 
signal-rich region with 55 < Mjj < 120 GeV/c 2 and 
two control regions with very little expected signal, 
Mjj < 55 GeV/c 2 and Mjj > 120GeV/c 2 . We com- 
pare the observed distributions to the predicted distri- 
butions in these three regions. In addition to the ma- 
trix element input quantities, we also check the mod- 
eling of the$r and mxiW) to validate the description 
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Figure 1: Mjj for events with EPD<0.25 (top left), 
0.25<EPD<0.5 (top right), 0.5<EPD<0.75 (bottom left), 
and EPD>0.75 (bottom right). 



of non-VF events. Finally we check ARjj, and PTjj to 
validate the modeling of the correlations between the 
two jets. All of these quantities are well-described by 
the simulation for our event selection. 

A small discrepancy in the modeling of Mjj itself 
can be seen when combining the two control regions; it 
is also visible in the low-EPD region of Fig. [TJ We as- 
sign a systematic uncertainty on the shape of the EPD 
in the W+jets background associated with this dis- 
crepancy by reweighting Monte Carlo events to make 
them agree with data. The weights are derived in the 
control regions and extrapolated through the signal 
region to avoid bias from the expected signal. This 
uncertainty has a negligible effect on the results, be- 
cause a large fraction of background events lie in the 
first few bins of the EPD distribution. Small changes 
in modeling of those background events do not signif- 
icantly change the shape of the EPD. 



torization and renormalization scale (Q 2 scale) used 
in the ALPGEN W/Z+ jets simulations, and the inte- 
grated luminosity. 



7. Results 

The observed and predicted EPDs are shown in Fig- 
ure We use a binned maximum likelihood fit of the 
observed EPD to a sum of templates. Systematic un- 
certainties are incorporated into the fit as nuisance 
parameters. We perform pseudo-experiments to cal- 
culate the probability (p-value) that the background- 
only discriminant fluctuates up to the observed re- 
sult (observed p-value) and up to the median ex- 
pected signal plus background result (expected p- 
value). We observe a p-value of 2.1 x 1CP 7 , corre- 
sponding to a signal significance of 5.4er, where 5.1a 
was expected. The observed WW + WZ cross section 
is 17.7±3.1(stat)±2.4(sys) pb. 



10 4 



CDF Run II Preliminary, L=2.7 fb 1 



^Data 
□ WW+WZ 
m W/Z+jets 
QCD 
Top 




Figure 2: Observed EPD distribution superimposed on ex- 
pected distribution from simulation. 



6. Systematic Uncertainties 

In addition to the uncertainty due to the Mjj mis- 
modeling described above, we consider several other 
sources of systematic uncertainty, taking into account 
their effect on both the signal acceptance and the 
shape of the background and signal templates. The 
uncertainty on the normalization of the backgrounds 
is considered part of the statistical uncertainty. The 
largest systematic uncertainty is due to the jet en- 
ergy scale (JES) uncertainty. Its effect is quantified 
by varying the jet energies in the simulated samples 
by ±lcr. We assign both an acceptance uncertainty 
and a shape uncertainty on the signal templates due 
to the JES. Other sources of systematic uncertainty 
considered are initial and final state radiation, parton 
distribution functions, jet energy resolution, the fac- 



8. Search Using Fit to 

A complimentary search for a resonance in the in- 
variant mass distribution has been carried out by fit- 
ting the Mjj shape to a sum of signal and background 
templates. The background should be smoothly 
falling while the signal should exhibit a peak at around 
80 GeV from the hadronic decay of a W or Z boson. 
This is a somewhat simpler search technique than us- 
ing matrix elements as described above, but is not 
expected to be as sensitive since it does not take ad- 
vantage of the full event kinematics. A larger data 
sample was used for this search, corresponding to 3.9 
fb -1 of integrated luminosity. 

The event selection for this search is somewhat 
different than what has been described for the ma- 
trix element analysis. The central lepton triggers 
are used, and offline events are required to contain 
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one electron candidate with Et > 20GeV or one 
muon candidate with px > 20GeV/c. We require 
£t > 25 GeV and at least two jets with Et >20 
GeV. The non-VF background is reduced by requir- 
ing i71t(W) >30 GeV/c 2 . Events are rejected if the 
two jets have A?j(Jl, J2) > 2.5. We also reject events 
whose hadronic vector boson candidate has pr < 40 
GeV; this selection improves the agreement between 
data and simulation and ensures a smoothly falling 
shape in the background dijet invariant mass in the 
expected signal region. 

On data, we estimate the signal fraction by per- 
forming a x 2 fit to the dijet invariant mass separately 
for the high pt muon and electron samples. Three 
Mjj template distributions are used in the fit. The 
first is the combination of W+jets, Z+jets, tt and 
single top production (referred to electroweak (EWK) 
backgrounds) where the relative normalizations are es- 
timated using predicted cross sections and efficiencies 
derived from simulation. The overall normalization of 
the EWK contribution is a free parameter in the signal 
extraction. The second template is the multi-jet QCD 
background template; its normalization is constrained 
by the fit to the-Er described earlier. The third tem- 
plate describes the signal and its normalization is a 
free parameter in the fit. In addition, the global nor- 
malization of the sum of our signal and background 
models is a free parameter. Figure [3] shows the data 
superimposed on the fitted templates after the elec- 
tron and muon samples are combined. 



CDF Run II Preliminary \L dt = 3.9 fb" 







~\ 






— Data 
^] EWK+QCD 
WW/WZ 
Fit 




























, , , 



50 100 150 200 

M„ GeV/c 2 



Figure 3: Dijet mass distribution in data superimposed on 
the predicted background and signal distributions. 

The systematic uncertainties in the fit to the Mjj 
distribution are treated in a similar way as in the ma- 
trix element search. In this case, however, the largest 
sources of systematic uncertainty are the shape uncer- 
tainties of the EWK and multi-jet background tem- 
plates. The jet energy scale uncertainty is also signifi- 
cant, as is the overall 6% uncertainty on the integrated 



luminosity. 

We estimate 428 ± 177 (stat) ± 56 (sys) events in 
the WW/WZ -> evjj sample and 650 ± 149 (stat) ± 
66 (sys) events in the WW/WZ — > nujj sample, giv- 
ing cross sections of 10.3 ± 4.2 (stat) ± 1.7 (sys) pb in 
the electron sample and 19.5 ± 4.7 (stat) ± 2.8(sys) 
pb in the muon sample. Combining the two decays, 
we estimate a total of 1079 ± 232 (stat) ± 86 (sys) 
WW/WZ — > Ivjj events, corresponding to a signifi- 
cance of 4.6 a and a cross section of o-yyw/wz+x = 
14.4 ± 3.1 (stat) ± 2.2 (sys) pb. This result agrees 
with the cross section measured with the matrix ele- 
ment technique as well as with NLO prediction. 



9. Conclusions 

In summary, we have observed WW + WZ pro- 
duction in the lepton plus two jets final state. We 
performed two searches: one building a discriminant 
using matrix element calculations and the other look- 
ing for a resonance on top of a smoothly falling di- 
jet mass distribution. The signal was observed with 
5.4ct significance with the matrix element search, and 
4.6ct evidence of the signal was found in the Mjj 
search. The WW + WZ cross section was mea- 
sured to be 17.7 ± 3.1(stat)±2.4(sys) pb and 14.4 ± 
3.1(stat)±2.2(sys) pb with the matrix element and 
Mjj searches respectively. These two results are in 
good agreement with each other and with the NLO 
prediction of 16.1 ± 0.9 pb. 
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